% Change file for DDTONE.WEB, computes a "sphere" (no, actually constant)
% Change file for DDTONE.WEB, changes to yield the Floyd-Steinberg method
@x
reset(bytes_in,input_name,'/B:8')
@y
@z
@x
@!t:eight_bits; {byte of input}
begin new_row[0]←0.0;
if i>mm then for j←1 to nn do new_row[j]←0.0
else for j←1 to nn do
	begin read(bytes_in,t); new_row[j]←(255.5-t)/256.0;
@y
@!x,@!y,@!z:real; {coordinates of input}
begin new_row[0]←0.0;
if i>mm then for j←1 to nn do new_row[j]←0.0
else for j←1 to nn do
	begin new_row[j]←0.5;
@z
@x
@<Choose pixel values and diffuse the errors in the buffer@>=
for k←0 to 63 do
	begin i←class_row[k]; j←class_col[k];
	while j≤nn do
		begin @<Decide the color of pixel |[i,j]| and the resulting |err|@>;
		for l←start[k] to start[k+1]-1 do
			begin u←i+del_i[l]; v←j+del_j[l];
			buffer[u,v]←buffer[u,v]+err*alpha[l];
			end;
		j←j+8;
		end;
	end
@y
@<Choose pixel values and diffuse the errors in the buffer@>=
for i←1 to 8 do for j←1 to nn do
	begin @<Decide the color of pixel |[i,j]| and the resulting |err|@>;
	buffer[i,j+1]←buffer[i,j+1]+err*alpha_e;
	buffer[i+1,j-1]←buffer[i+1,j-1]+err*alpha_sw;
	buffer[i+1,j]←buffer[i+1,j]+err*alpha_s;
	buffer[i+1,j+1]←buffer[i+1,j+1]+err*alpha_se;
	end
@ @<Const...@>=
alpha_sw=0.1875; {3/16}
alpha_s=0.3125; {5/16}
alpha_se=0.0625; {1/16}
alpha_e=0.4375; {7/16}
@z
%alpha_sw=0.166666667; {3/18}
%alpha_s=0.277777778; {5/18}
%alpha_se=0.055555556; {1/18}
%alpha_e=0.388888889; {7/18}